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Abstract 



CP violating effects in top quark pair production at the future 2 TeV 
pp Tevatron and 14 TeV pp LHC colliders are investigated. We study 
three kinds of CP violating sources: the supersymmetric CP-odd phase 
of the stop trilinear soft breaking term: arg{At), the CP-odd parameter 
in two-Higgs doublet extensions of the standard model(2HDM), and the 
model-independent top quark chromoelectric dipole moment(CEDM), 
respectively. Optimal observables as well as simple observables are 
used. We find that it is possible to observe CP violating effects from 
arg{At) in top quark pair production at the 2 TeV Tevatron with ~ 
30/6~^ integrated luminosity when nig ~ 200 GeV. If the experimental 
systematic errors are sufficient small, the LHC with ~ 150/6~^ can put 
a limit of order 10^^ on the phase arg{At) and the CP-odd parameter 
in 2HDM by using optimal observables. The CEDM of the top quark 
can be measured to an accuracy of 10~^^ cm Qs at the Tevatron and 
few X 10-20 cm at the LHC. 
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I. INTRODUCTION 



At the future hadron colliders such as the upgraded Tevatron with y/s = 2 TeV 
and the CERN LHC with -^i = 14 TeV, the annual top quark pair yields are about 
7 X 10^ (with integrated luminosity /£ = lO/fo-^) and 7 x 10^ (with JC = 100/6"^), 
respectively. The large numbers of top quark pairs allow us to do precise measurements 
on physical quantities associated with the top quark. Among them are the production 
cross sections, the top quark mass and the discrete symmetry properties. Due to 
the QCD uncertainties and the experimental systematic errors, the precision of cross 
section measurement is only about 5% — 6% The discrete symmetry properties 
such as parity nonconservation and CP violation do not suffer from QCD background 
uncertainty and their accuracies of measurements depends mainly on the statistical 
errors provided the experimental systematic errors are sufficient small. Therefore the 
discrete symmetry properties can be measured more precisely than the cross sections. 
New physics which has no observable contributions to cross section may have observable 
effects in parity or CP violation. Since the Standard Model(SM) contributions to parity 
violation and CP violation are small, possible observed large effect of them will reveal 
new physics. 

In a previous work 0, we have shown that the minimal supersymmetric(SUSY) 
extension of the standard model (MSSM) |^ gives observable parity violating effects at 
the Tevatron while its corrections to the production cross section are within the QCD 
theorectical uncertainties. In this paper, we shall concentrate on the CP violating 
effects induced by non-standard model interactions. We investigate three kinds of 
CP violating sources: the supersymmetric CP-odd phase of the stop trilinear soft 
breaking teTm:arg{At), the CP-odd parameter in two-Higgs doublet extensions of the 
standard model(2HDM), and the model- independent top quark chromoelectric dipole 
moment (CEDM), respectively. In the MSSM, CP violation exists in strong interaction 
and in the 2HDM, CP violation can have strong Yukawa couplings due to the heavy 
top quark mass. Therefore, both can produce possible large effects. It is also useful 
to study CP violation in a model-independent way when we do not know what is the 
new physics. CP violation in top quark pair production at hadron colliders is also 
studied in Refs. [|]- In Ref. 0, SUSY QCD CP violating effects are studied in 



the gg ti process by using a charge energy asymmetry observable which is only 
sensitive to the imaginary part of the loop integrals. In this work, we use the optimal 
observables as well as naive observables constructed from the final state momenta. 
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Possible large CP violating effects in 2HDM and the methods of observing them in 
top quark pair production at hadron colliders are studied in Refs. [|T^. We extend 



those studies by applying optimal observables. In Ref. 0, the method of extracting 
real top quark CEDM in the reaction gg — > ti is studied. It is found that the optimal 
observables are particularly effective. We include here an imaginary part of CEDM and 
the reaction qq —>■ ti at the Tevatron. Furthermore, we use the exact amplitudes of 
QdilQ) —^ti^ hli Uij}l2 ( bqiq[bq2q2)- In Ref. P], the top quark spin in its rest frame 
is taken to be in the direction of the lepton. This is a kind of approximation, althogh 
the lepton is a good analyzer of the top quark spin, because in the top rest frame, the 
lepton momentum has the angular distribution proportional to 1 + cos ip with ip being 
the angle between the top spin and the lepton momentum. The possibility of using 
polarized proton was studied in Ref. . 



In Sec II, we describe the models and the calculations. The methods of extracting 
CP violating effects are given in Sec. III. In Sec IV, we present our results, discussions 
and conclusions. 

II. MODELS AND CALCULATIONS 



A. CP violation in MSSM 

In the MSSM, two possibilities can induce CP violation in top quark interactions: 
the complex phase in Higgs mass parameter fi, and the complex phase in scalar top 
supersymmetric soft breaking trilinear coupling At. The experimental limit on the 
Neutron Electric Diploe Moment(NEDM), dn < 1.1 x 10~^^e — cm [|15|, places a severe 
constraint on the phase of /i. Therefore, the only significant SUSY CP-odd phase in 
associate with the top quark is arg{At). In Ref. [|16|, it is argued that the phase arg{At) 



is not strongly constrained by current experiments and the effects in single top quark 
production and decay are studied. In this work, we shall assume arg{fi) = and let 
arg{At) to be a free parameter of no a priori constraints. 

The parameter arg{At) enters in the scalar top quark mixing. The mass eigenstates 
ti and ^2 of scalar top quark are related to the current eigenstates ti and t/j by 

ii = iL cos 6t + in. sin Ote''^', t2 = sin 6*46*^* + in cos Ot (1) 

The mixing angle 9t, phase Pt as well as the masses m^^ ^ can be calculated by 
diagonalizing the following mass matrix M 
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(Ml m,ml^ \ 
* \mmLR Ml I ' 
1 2 

Ml = ml + + (2 - 3 si^' ^w) cos(2/?)m| , 
2 

= ml +m^ + - sin^ 6^1^ cos (2/5) m| , 
"^Lij = -/i cot P - At , (2) 

where m?^ , m?^ are the soft SUSY-breaking mass terms of left- and right-handed stops, 
tan P = V2/V1 is the ratio of the vacuum expectation values of the two Higgs doublets. 
From Eqs.(l)(2), we can get the expressions for rnf^^, 9t and (3t : 

1 



mf = - 



M4+M4Ty'(M^-M^)2 + 4m,>iR|2 , (3) 



i^^9t = — !\ , l\ ^cosA , (4) 

mt(—/i cot p — \At\ cost) At) 



t^^A = ' I , , , (5) 





At 


sin 6*^^ 


cot /3 -|- \ At 


cos 6 At 



where 6 At = 0,1^ 9 (At) . 

In the presence of squark mixing,the strong squark-quark-gluino interaction La- 
grangian is given by 



Lgqq = - Q sT^k^k^ai - 6l75)gij + (^2 - &275)<?2j]5'a + H.C. , (6) 

where gs is the strong coupling constant , are SU{3)c generators and for top quark 
fli, bi, 02, ^2 are given by 

ai = ^ (cos 6^^ - sin 9te'^' ) , h = — ^ (cos 6't + sin Ote''^' ) , 
v 2 V 2 

02 = — ^(cos^t + sin^^e-*^'), 62 = — ^(cos^f - sin ^^e"*'^*). (7) 
V 2 v 2 

The above interactions enter in the virtual corrections to the main production pro- 
cesses of ti at hadron colliders: qq —>■ ti and gg —>■ ti. There are also weak squark-quark- 
neutralino and squark- quark- chargino interactions. Since their coupling constants are 
an order of magnitude smaller than the strong SUSY QCD squark-quark-gluino inter- 
action, we shall not consider them here. However, they may be important when the 



main processes are weak interactions |]I6 



In Fig.l(a)-(e), the Feynman diagrams of the QCD tree level and SUSY QCD 
virtual corrections of the process qq —>■ ti are given. The corresponding Feynman 
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diagrams of gg ti are presented in Fig.2(a)-(m) (the u-channels of Fig.2 (b),(f)- 
(m) are not depicted). The dashed hnes in the loop stand for scalar quarks, while 
the solid lines for gluinos. In Ref. [0, CP violating effects are studied using a charge 
energy asymmetry observable which is only sensitive to the imaginary part of the loop 
integrals. Therefore, there are no contributions to the charge energy asymmetry from 
Fig.2 (f)-(j),(m). 

The one loop scattering amplitudes of qq ti and gg ti were already presented 
in Refs. [p!7Hl^ for calculating the total production rates of ti pairs. To calculate the 
CP violating effects in the ti system, additional renormalized amplitudes are needed. 
In terms of the tree-level amplitude, Mq , and the next-to-leading order SUSY QCD 
corrections, SM"-, the renormalized amplitudes of ad ti {a = q, g) at the one-loop 
level may be written as 

M" = Mq" + SM". (8) 

where SM'^ can be decomposed into two parts: 6M"'^ which contains even combination 
of 75 and e^,^po-, and 6M"-^ containing odd combination of 75 and e^^^po-. The symmetry 
breaking effects are contained in ^M""^ which has no contributions to the total cross 
sections at next-to-leading order, while 6M"'^ will contribute to the total cross sections. 
We shall assume that SM""^ is small enough to be within the 5% — 6% uncertainty 
therefore is neglected in our calculations. We also discard terms in SM""^ which give 
only parity asymmetry 

Let us denote the momenta of the initial and the final state particles as 0(^4)0(^3) —* 
ti{P2)ij{pi)- We may use, as a further short-hand, the notation for a = g that Ui = u{pi) 
[vi = v{pi)) denotes the Dirac four-spinor corresponding to the momentum and spin 
of particle (anti-particle). When a = g we use = e{pi) for the gluon polarization 
function. In this notation, the tree level amplitude for a = q can be written as 

MS = igl{T^;r^Jv,^^u,u^^,v,/s, (9) 

where s is the invariant mass of the ti pair. 

The tree level amplitude for a = g is composed of three different production 
channels(s-,t-,u-channel) as following: 

M'o' = -tg!itfabcn,MP2)Mpi)/S 

= -igliT'^T^ - T^T'')jiM'^, (10) 
MS' = -igl{T'T'^)^iu{p2)iM + m,)i,v{p^)/{i - ml) 
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= -igl{T'T'^)^,Ml (11) 

= -tgl{T'^T%,M^, (12) 

where q = P2 — Pa, T'^ is given in the Appendix A. 

To calculate the CP violating effects induced by the SUSY QCD effects, we follow 
the method presented in Ref. [|5[ [0, in which the amplitudes were calculated nu- 
merically using the helicity amplitude method. To obtain the renormalized scattering 
amplitudes, we adopt the dimensional regularization scheme to regulate the ultraviolet 
divergences and the on-mass-shell renormalization scheme to subtract the divergences 



The SUSY QCD corrections to the scattering amplitude SM"-^ for a = q arise 
from the vertex diagram, the box diagram, as well as the crossed-box diagram. The 
renormalized amplitudes can be written as 

(5M''^ = SM"" + 6M^^ + 6M^^, (13) 

where SM'^^ is the vertex corrections of Fig. 1 (b)(c), and SM^^ and SM^^ are the con- 
tributions from the box diagram and crossed-box diagram of Fig.l (d)(e), respectively. 
The results for these separate contributions are, 

6M^^ = igliT^.TtJHh - ^i)uM'D'lMs, (14) 

+U2U4V3:^475Vi + M275M4t^3|^4t^l] , (15) 



5M^^ = ^<5M^^(p3 ^ Pa). (16) 

The form factor D** corresponding to the top quark CEDM and the form factor f^^ 
are given in the Appendix A. 

The SUSY QCD corrections to the scattering amplitude 5M^^ of gg — > tt can be 
written as 

SM^f = -igl[5M+0^^^ + 5M-0]7) + 5M\5abU , (17) 

with 
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3 

-(5M^<'^^'* + 5M^°^^'") (19) 

In the above expressions, the superscript t,u stand for t-,u-channeL 5M^^ is the 
s-channel vertex corrections from Fig. 2 (c)(d), 5M^^ from (e), 5M^^^^'* the self energy 
correction from (f), 5M^^'* from (g)(h),5M^2,t (j^^j)^ ^m^^^^'* from (k), (5M^°^2,t 
from (1) and 5M^°^^'^ from (m). In the following we only give the explicit results of 
the s-channel (no crossed diagram) and t-channel contributions. The u-channel results 
can be obtained by the following substitutions: 

Ps ^ Pi, a ^ b, i ^ it. (20) 

All S-, t-channel terms in SM~, SM^ can be written as following according 

to their Lorentz structures: 

= e^;^M2 [fiQi,. + fflf^lu + fj^'pi^Plu + fsPli.P2. + fgP2f.Plu (21) 
+fwP2fiP2u + fl2,i>APlvl„ + fu^muli^ + &4P2ulix + fie,hP2nlv]lbVl , 

where X = si, s2, self, vl, v2, boxl, box2, box3, respectively. The 10 form factors 
corresponding to each diagram are given in the Appendix A. They are not all indepen- 
dent when we sum over all possible channels. CP-odd property of 6M^ requires 

fn — fw, fu — /i5 ' fr^ + = /{^ . (22) 

Those relations arc verified by our explicity formula of the form factors. We find that 
all form factors are proportional to Xcp = 2Im{aib{) = sin 2^^f sin /3f . 

The color sum of the amplitude square including the next-to-leading order correc- 
tion is: 

^{\M^\'' + 2Re(M^5M'^^)} (23) 



color 

9t{lM\'+mo\' 



+^Re{M+dM+^) + 6Re{M^SM-^) + 8Re{M+dM^^)} , 
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where 



Mo+ = M* + Mo", Mo" = M* - Mo" - 2Mo", (24) 

B. CP violation in 2HDM 

In ordinary 2HDM, there are three neutral physical Higgs bosons, namely two CP- 
even scalars H, h, one CP-odd pseudoscalar A. CP violation in the scalar potential 



22| induces mixing of the CP-even and CP-odd neutral Higgs bosons, thus leading to 
three physical mass eigenstates \(f)j > {j = 1,2,3) with no definite CP parity. Their 



Yukawa couplings to the top quark can be written as(in the notation of |]23 



Ly = -m<(V2G^)i/2 ^ ^ m,t75)% (25) 

where Gp is the Fermi constant, and 

o-jt = d2j/ sin/? , djt = —d^j cot P , (26) 
, d^j are the matrix elements of a 3 x 3 orthogonal matrix which describes the mixing 



of the neutral states p3|. We assume that the two heavier Higgs bosons(j = 2, 3) may 
be neglected and define a = au, d = au, (p = (pi. The strength of CP violation is 
proportional to 2/m[aif(— iau)*] = 2ad = —2'ycp{lcp is defined in [0). 

The one-loop Feynman diagrams of Higgs boson contributions to the processes qq 
ti and gg tt can be represented by Fig. 1(b) and Fig.2(c)(f)(g)(i)(k)(n), respectively. 
Now the dashed lines in the loop stand for Higgs bosons, while the solid lines for top 
quarks. 

The one-loop amplitudes can be easily obtained from those of the MSSM by the 
method given in the Appendix A with no changes on 5M'^'" , 5M^ , 5M~ and the 
modification on 5M^ which comes from Fig. 2 (n): 

5M' = etelu2{a\n^g,, + n'p,,p,.) - fl'ie^.p^plvli, (27) 
+ad[i{fl''g^,y + /|>3/.P4i/)75 - fl'' (^^lupcJPlpV\}vl/ (s - mj + im^V^) 

where is the Higgs boson width and 3 given in the Appendix A. 5M^ is 
composed of two parts: 5M^^, which is CP-even and contains terms of a?, d?; 5M^° , the 
CP-odd term proportional to ad. We keep the CP-even term because when > 2mt, 
it may be important. The ^Rei^M^ 5M^^) term in the color sum of ( ^31) should be 
replaced by 
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8Re{M+5M^^) + 2A{\5M^^\'^ + \5M^°\'^) + A8Re{5M^^5M^°^) . (28) 

We compute by using the formula and parameters of Ref. |12| 

C. CP violation from model-independent top quark CEDM 

When the particles in the loops are heavy compared to the external particle mo- 
menta, it is convenient to describe the loop induced interactions by effective interaction 
Lagrangian. Even though when the loop masses are not too large, a model-independent 
study can give us knowledge about the sensitivities of given colliders. That approach 
is particularly useful when we do not know the underlying new physics. We assume 
the following additional top-quark-gluon effective interaction: 

Ld = -^ta>^^^,F^^,TH , (29) 

where a^" = (i/2)[7'^, 7^], Fj^^ is the the gluon field strength. The coefficient df is the 
top quark chromoelectric diploe moment (CEDM) which we assume to have imaginary 
part as well as real part. We denote it as + idf. 

The CEDM dt contributes to the CP violating amplitudes through the diagrams of 
Fig. 3. Let us define dt = dt/ Qs = d^ + idf. The contribution of Fig. 3 (a) to qq —>■ ti can 
be obtained by simply replacing with idt in Eq. ([l^) and setting 6M^^ = 6M^'^ = 0. 
We denote the contributions to gg ti of Fig.3(b)-(e) as 5M^ 5M^ ^M"-*, 5M'''"(for 
the crossed diagrams), respectively. Then we have 

5M+ = dM"^^ + (5M^'" , 

5M- = -26M' - 26 M'' + 5M^'* - 5M^'" , 

6M^ = . (30) 

SM'^, 6M^, can also be written as the form of Eq. (|2l|) with the constant form 

factors given in the Appendix A. 

III. EXTRACTING THE CP VIOLATING EFFECTS 

Since SM"""^ contains odd combination of 75 and e^^pa, its interference with the 
tree level amplitude will be zero if we sum up all the initial and final state spins or 



The misprint of Tz is corrected. 
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polarizations. To see the CP violating effects, the spins of at least some of the particles 
must be identified. Because the initial state spins are difficult to be determined, the 
helicities of the final ti must be inferred from their decays. If we assume the SM 
couplings of the top quark to its decay products and set = 0, all decay products 
of the top quark will have left-handed helicities. That means the spin information of 
the top quark can only be transfered to the momentum correlations among the decay 
products. 

We assume the SM decay of the top quarkQ and apply the narrow width approxi- 
mations of the top quark and W-boson propagators: 

1 TT 

—S{q^-ml). (31) 

where Y stands for top quark and W-boson, Fy is the width of Y 
The parton level cross section for reaction ad ^ ti ^ blf ui-Jjl^^ 
be written as: 



da., = (32) 



where 7 = y 1 — Amq / s and 

^ ,j _ + - « (rn? - r,4)V< . (33) 

dQ'^+{dQ'^^-) is the solid angle element of Vr^(iy~) in the rest frame of the (anti) 

top quark, dQ'+{dQ' ) denotes the solid angle element of /j^(/^) in the rest frame of 
2 

W^{W~), \M\^ is the average amplitude square excluding the top quark and W-boson 
propagators after the decays of the top quarks: 

W\l = S-' E {|Mo"P + 2i?e(MoW^t)}^ (34) 

color, spin 

where l/Sa is the color, spin average factor: Sg = 36 and Sg = 256. In our calculations, 

2 — 

|M|^ is easily obtained from the amplitude of ad tt by the following substitutions: 
U2 -^«67m(1 - 75) (1^2 + mt)u^^^^{\ - 75)wi, , (35) 



The top quark decay asymmetry in the MSSM is found to be small |16l [24| 
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where g is the weak SU{2) couphng constant. The above expresssions are calculated 
numerically. 

The hadronic cross section is obtained by convoluting (|32D with parton distribution 
functions: 

d(Tpp(p)= I dXi f dX2[f^iXi)fi^^\x2) + (p Pip))]daaa/{1 + Sa,a) , (36) 

where /f^^^ are parton distribution functions(PDFs). 

It has been studied extensively in the literature how to pick out the momentum 
correlation information contained in 2Re{M^5M'^^'^) and therefore the CP violating 
effects by looking at the expectation values of the CP-odd observables constructed 



from the momenta of the top quarks and their decay products |0- [|T6| ||2^- 

In this study, we adopt the following simple observables which are constructed from 

observed momenta and can be easily used by experimentists: 

= Ei^ - Ei^ , (37) 

M=Pt- Pi+ - Pt ■ Pi~ = Ot - Qt , (38) 

T2 = {Pb - Pb) ■ {Pi+ ^Pi-) , (39) 

r ^jiuapPl + Pl-PbPl /.^N 

{pi+ ■pi-Pb-Pir'^ 

= —r;^/-iPi+xpi-)P-iPi^ -Pi~) ' (41) 

mi\Fy 

where all momenta are in the laboratory frame, Ei+{Ei-) is the energy of = 

h = h = e, fi, here we do not distinguish e, /i), the subscripts of the momenta 
denote the corresponding particles, P is the momentum of the proton in the case of 



pp collision. Ai, A2, and T2 are studied in Ref. [12|. /2 and Ol are used in Refs. 



|Tl| , respectively. Because Ai, A2 are T-(time reflection |]29|) even, they are 
only sensitive to the absorptive parts (imaginary parts) of the loop calculations and 
df. On the contrary, T2, /2, and Ol are T-odd and only sensitive to the dispersive 
parts(real parts). Ai, T2, f2 and Ol all require the events that two top quarks decay 
semileptonically. A2 uses the events that one top quark decays semileptonically and 
the other hadronically. Among the above observables, only /2 is Lorentz invariant. 

While all the above naive observables use only parts of the information, the opti- 
mal observables studied in Refs. 0]- use the full information in 2Re{MQSM'^^'^). 
Therefore they are the most effective ones. 

In the case of model independent top quark CEDM, 2Re{MQSM"'"^^) contains two 
terms which are proportional to d^ and d(, respectively. In the MSSM , it depends on 
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the particle masses in the loop as well as a multiplicative constant Xcp- However, since 
our main goal is to search for CP violation induced by this Xcp-, we must assume all 
the masses in the loop are known. In the 2HDM, from the first term in Eq.(p8D which 
is CP-odd, one can separate a factor '-jcp = —ad. Although in the resonant region, 
contributions of the other three terms ( belonging to higher than next-to-leading order) 
may be large, they are still overwhelmed by the tree level and the next-to-leading 
order contributions. As an approximation, we shall drop them in the definitions of 
the optimal observables. Therefore we can always separate a constant (denote it as 
A = Xcp, Icp, df, d() from 2Re{MQ6M"-^'^) in all the models considered. 

Apart from some common factors in (|36|) , the hadronic cross section can be written 

as 



d^PPiv)=llj E {|Mol2 + A2/?e(MoW^t)} (42) 

O'—QtB color, spin 



where d^ denotes the phase space. In the following, 2Re{MQ6M"-^^) is calculated by 
setting Xcp = 1, 7cp = ^, dt = ^ and dt = i in the form factors in each model, 
respectively. The optimal observable is defined as 

E E 2ReiMSSM''^^)[fP{xi)jf''\x2) + {p ^ pip))] 

Q a— QiS color, spin f43) 

E E \MS\VSix,)fi^^\x2) + ip^ Pip))] 

O'—QtS color, spin 

The above observable depends on the parton distribution functions. It is inconvenient 
for practical use. In the pp collision, it is not a CP-odd observable due to the asym- 
metry of quark PDFs (cf. Appendix C). Because at the 2 TeV Tevatron, the main ti 
production mechanism is qq tt, and at the 14 TeV LHC, the main process is gg ti, 
we can neglect one process at each collider. We consider only qq —>■ tt at the Tevatron 
and gg tt at the LHC in all the following calculations. Then we have the optimal 
observable 

E 2i?e(Mo"5M"^t) 



E |Mo"|2 ^ ' 



color, spin 

color, spin 

Since the neglected process consists of only about 10% of the total cross section, Eq. 
(|4^) deviates from the truly optimal observable by at most 20%. Oi has the property 
of Lorentz invariance, one can calculate it in any frame. It has the same symmetry 
property as 2ReiM^6M''^^) which is CP-odd. In the MSSM and 2HDM, d has no 
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definite T parity and depends on the loop particle masses. We note that the optimal 
observable defined for is only sensitive to independent of d^. The same holds 
true for d^. This is because the two terms proportional to (if and df have different 
discrete symmetry T. 

To calculate Oi, we need to know all the momenta of the initial and final state 
particles. That can not always be achieved. It is still useful because it can provide us 
with the upper limit signal to noise ratios of any other CP-odd observables. Let us 
first look at the case that two top quarks decay semileptonically. The two neutrino 
momenta p^-^, p^^ will be unknown. But they may be determined indirectly by the 
following eight equations [0: 

pI, = PI2 = 0, {Pu^ + Pi+f = {Pu2 + Pi-f = ml,, (45) 
{Pui + Pi+ +Pbf = ml, + pi- +pif = ml, 

{Pu2 ~^ Pl~ Pi)) transverse iP^i ~^ Pl^ ~\~ Pb)transverse ■ 

Similar situation exists in the study of rf production in e~^e~ collision [Q. The 
above equations give rise up to fourfold solutions(see Appendix B). Considering this 
ambiguity, we can define a modified optimal observable: 

E 2Re{Mg5M''^^)r], 

^ color, spin,i / a n\ 

°' = E |A/S|=,, ■ '"'^ 

color ,spin,i 

where the sum i is over all possible solutions of the neutrino momenta. 77^ = 
Ti/'^J/a (^^D/f ^^^^2) + (p ^ P{p))\ comes from the tt phase space, fiux factor and 



PDFs due to different momentum reconstruction(cf. Eqs. (p2D and (|36D ). There may be 
possibility that the reconstructed initial parton energy exceeds the proton(antiproton) 
energy. That kind of reconstruction should be discarded in the calculations. 0/ also 
depends on PDFs. For practical use, we define a non-optimal observable as an approx- 
imation: 

0[ = Y.Oi. (47) 

i 

Again the sum i is over all possible solutions of the neutrino momentum^. 



^One can also define 0[ by setting r/j = 1 in (^), the numerical difference between the two 
definitions is minor 
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We now consider that one top quark decays semileptonically and the other hadron- 
ically. The missing neutrino momentum can be fully reconstructed [Q. But because 
we can not distinguish quark and antiquark jet, we still have ambiguity of twofold 
uncertainty. When two quarks decay all hadronically, there is fourfold uncertainty. We 
define therefore alternatvely the optimal observables : 

E 2Re{Mg6M''^^) 

^ color, spin, j /Ar^\ 

°' = ^liw — ■ 

color,spin,j 



color , spin, j' 



= — ^-'^ , (49) 

color, spin, j' 



where the sum j in O2 is over the two possible assignments of the jet momenta, and 
the sum j' in O4 is over the four assignments 

The statistical significance of an observable O can be described by the signal to 
noise ratio r defined as 



r=<0> /V< 02 > , (50) 
where < O >,< > are the expectation value of O, O^, respectively: 

< >= Z^!^ . (51) 

It is interesting to note that for the optimal observables and unit A, we always have 
< O >=< >. Care must be taken in calculating the r of yl2- Because A2 is the 
difference of two observables (Of and Oj) which are calculated using different events(i.e. 
different independent distribution functions), we have < >=< > + < >~ 
2 < 0| >. If the experimental error comes only from statistics, the number of events 
Nevent needed to observe CP violating effects at la level (68% C.L.) satisfys \r\ > 
N event OT Nevent > l/r^. To rcducc the Statistical errors, one can combine the 
measured results of the three decay modesdeptonic-leptonic, leptonic-hadronic and 
hadronic-hadronic modes |T5[. Assuming their corresponding number of events are 
Nii,Nji and Njj, respectively, then we can define a combined signal to noise ratio: 



Vc = ^{Nurl + N^irl + N,fl)/N, (52) 

where ri, r2 and are the signal to noise ratios of the observables which use leptonic- 
leptonic, leptonic-hadronic and hadronic-hadronic events, respectively. = Nii + Nji + 
Njj is the total number of events of the three modes. Note that depends only on the 
ratios Nu/N, Nji/N and Njj/N, not on A^. The signal is detectable at la level when 
\rc\ > 1/VN 
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IV. NUMERICAL RESULTS AND CONCLUSIONS 



We first check our calculations with QCD gauge invariance in the process gg — > ti. 
That can be done by replacing ea, 64 with ps, p^. We find that the correction amplitude 
is consistent with gauge invariance. Then we check our calculations with known results 
in the literature. Our parton level results are all in agreement with Refs. [0] |T2|. By 



using ^/s = 40 TeV, rrit = 160 GeV and 10'' sample of leptonic events, we find the Oi 
limit on df is 2.8 x 10~^° cm gg which is very close to Xmin in Ref. P]. That means that 
taking the top quark spin in its rest frame to be in the direction of lepton momentum is 
feasible. We can also reproduce the results of O^, in Ref. 0] with /i^ = by including 
both qq ti and gg ti processes in the calculations. 

As mentioned previously, in the following calculations, we consider only qq ti 
at the Tevatron and gg —>■ ti at the LHC The parton distribution functions of MRS 
set A' |35| with scale = are used^. rrit is take to be 176 GeV. To look at the 



largest possible effects, we set the CP- violating parameters Xcp and •ycp to be of order 
1, namely, Xcp = l(this needs Of = tc/A, cf. Appendix A), jcp = 1- In the MSSM and 
2HDM, we treat the SUSY particle masses and Higgs boson mass as free parameters 
allowed by current experiments. We assume that all the squarks except for the light 
stop ti to be degenerate. The light stop mass is required to be above 50 GeV ^7^. We 
choose gluino mass rrig and Higgs boson mass > 100 GeV. 

The signal to noise ratio results at the 2 TeV pp Tevatron and at the 14 TeV pp 
LHC are summarized in Table I. -Table V. We do not present the results for 2HDM at 
the Tevatron because the effects are too small unless the Higgs boson mass is less than 
100 GeV. We denote the signal to noise ratio r of an observable O as r(0). All the 
tables show that r{0'i) is slightly smaller than r(Oi) and r{02) is about 3/4 of r{Oi), 
r(04) about 1/2 of r(Oi). Therefore, 0[ is a good approximation for Oi. 

Table I. shows the results at the Tevatron from five sets of SUSY parameters. We 
see that the naive observables have signal to noise ratios all < 1%. It is difficult to 
observe such small effects at the Tevatron. The optimal observables, on the other 
hand, have r > 1% as long as the gluino mass is around 200 GeV Q. They are about 



■^It should be noted that all the results are insensitive to the PDF and choices. CTEQ4M 
gives similar results. 



^The value is close to mt, so that the gluino threshold is close to the top quark one. The 
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3 — 10 times more effective than the naive ones. The combined signal to noise ratios 
of 0[fi2 and O4 are all between these of O2 and O4. Because A*";; is small compared 
with Nji and Njj, the combined r is just a weighted average of these of O2 and O4. 
Assuming we can obtain 30 fh~^ integrated luminosity, then the total number of fully 
reconstructed tt events is about 3.5 x 10^ ||1|]. It is reasonable to assume that we can 
have purely hadronic events Njj = 2 x 10^, hadronic-leptonic events Nji = 1.2 x 10^ 
and purely leptonic events A^^^; = 0.2 x 10^. The corresponding 1 a level statistical 
errors are rjj = 0.71 x 10'^, rji = 0.91 x 10~^, ru = 2.2 x 10~^. The combined error is 

0. 54%. In the Table, O'l has r > ru only when m^-^ is around 50 GeV and rrig is around 
200 GeV. We always have r(02) > Vji, r(04) > rjj when rrig ~ 200 GeV. Therefore, 
it is possible to detect the CP violating effects by using optimal observables O2 and 

04 when rrig is around 200 GeV. The effects are detectable by combined measurement 
when rrig is in the range 100-300 GeV. 

Table II. is the results at the LHC in MSSM. It shows the same features as Table 

1. However, at the LHC, the numbers of events are much larger than those at the 
Tevatron. With 150 fb^^ integrated luminosity, we can assume Njj = 10'', Nji = 
6 X 10^ and A^^; = 10^. We further assume the experimental systematic errors are 
below the statistical ones. In addition, there are also theoretical uncertainties coming 
from possible non-CP violating contaminations at the pp LHC. Because the initial pp 
state is not a CP eigenstate, CP conserving interactions can produce CP asymmetry 
effects in tt final state. We present a general analysis of the contaminations to the 
Lorentz invariant observables in Appendix C. We find that within the framework of 
parton model, there are no contaminations to these observables. Ai, A2 and T2 are 
discussed in Refs. [|l^. They are well below the signals. Therefore, the 1 a level 
errors are rjj = 0.32 x 10^^, rji = 0.41 x 10"'^, ru = 1.0 x 10^^. The combined error is 
0.24 X 10^^. The naive observable T2 is on the margin to be detectable. Ai and /2 are 
better than T2. They can be used to observe Xcp at few xlO~^ when rrig is within the 
range 100 — 400 GeV. The observable O'l is 2-5 times better. All r{02) > 10 rji and 
^^(04) > 10 rjj for rrig ~ 100 — 300 GeV. We can pin down Xqp to 10^^ by using these 
optimal observables. 

In Table HI., we give the results of signal to noise ratios in 2HDM at the LHC. It 
is obvious that Ai and T2 are only detectable and can not be used to put limit on 7cp. 
/2 is 2 times better which may be used to limit 'Jcp to (3 — 4) x 10^^. All the optimal 



CP violating effects are large due to the threshold effects. 
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observables have r large than about 10 times statistical errors. Therefore, they will 
put a limit of order 10~^ on ^cp- 

Table IV. and V. are the results of model- independent top quark CEDM. An overall 
feature of the two tables is that both at Tevatron and LHC, the accuracies of dl are 
better than those of df. With the above assumed numbers of events at the Tevatron, 
the best limits on df and dl are 2.4 x 10~^^cm Qs and 1.1 x 10~^^cm gig, respectively, by 
using 02- Also by using O2 at the LHC, we can obtain the hmits of 5.2 x 10~^° cm Qs 
and 2.5 x 10~^°cm Qs on df and dl, respectively. 

To summarize, we have studied CP violating effects in top quark pair production at 
the future 2 TeV pp Tevatron and 14 TeV pp LHC colliders. Three kinds of CP violating 
sources:the SUSY CP-odd phase of the stop trilinear soft breaking term, arg{At), 
the CP-odd parameter in 2HDM, and the model-independent top quark CEDM are 
investigated. Optimal observables as well as simple observables are used. The optimal 
obervables are usually 2 — 10 times more effective than the naive ones. It is possible 
to observe CP violating effects from arg{At) in top quark pair production at the 2 
TeV Tevatron with ~ 30/6"^ integrated luminosity when rrig ~ 200 GeV. If combined 
measurement is applied, the range of gluino mass in which the CP violating effects are 
detectable is ~ 100-300 GeV. The LHC with 150/6"^ can put a hmit of order 10"^ on 
\cp (therefore on arg{At)) in MSSM and 7cp in 2HDM by using optimal observables 
provided the experimental errors are sufficient small. The CEDM of the top quark can 
be measured to an accuracy of 10~^^ cm Qs at the Tevatron and few xlO~^° cm Qg at 
the LHC. More accurate measurement on dl can be obtained than on d^ with given 
number of events. 
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TABLES 



TABLE L 

Signal to noise ratio r in pp — > + X at the 2 TeV Tevatron in the MSSM with 
\cp — 1, for five sets of SUSY parameters labeled by {mi^,mi^ — mq,mg) GeV. The 
combined results are for O2 and O4. 





Ai 


A2 


T2 


h 


Oi 


o[ 


O2 


O4 


combined 


(100,500,100) 


0.25% 


-0.23% 


-0.16% 


-0.16% 


0.82% 


0.79% 


0.62% 


0.45% 


0.54% 


(100,500,200) 


0.59% 


-0.42% 


0.08% 


0.12% 


1.82% 


1.73% 


1.41% 


1.04% 


1.18% 


(100,500,300) 


0.12% 


-0.17% 


0.13% 


0.11% 


0.78% 


0.74% 


0.56% 


0.40% 


0.49% 


(50,500,200) 


0.80% 


-0.58% 


0.08% 


0.14% 


2.49% 


2.37% 


1.93% 


1.43% 


1.69% 


(100,1000,200) 


0.58% 


-0.42% 


0.11% 


0.15% 


1.77% 


1.68% 


1.36% 


1.01% 


1.19% 



TABLE IL 

Signal to noise ratio r m pp ^ ti + X ai the LHC in the MSSM with \cp = 1, for six 
sets of SUSY parameters labeled by {mf^,m^^,mg) GeV. The combined results are for 
O'l, O2 and O4. 





-4i 


T2 


12 


Ch 


0/' 


O2 


Oi 


coiiil)ined 


(100,500,100) 


-0.12% 


0.11% 


0.22% 


0.92% 


0.81% 


0.71% 


0.53% 


0.62% 


(100,500,200) 


0.41% 


-0.13% 


-0.35% 


1.82% 


1.61% 


1.41% 


1.05% 


1.23% 


(100,500,300) 


0.17% 


-0.13% 


-0.29% 


0.85% 


0.72% 


0.62% 


0.44% 


0.53% 


(100,500,400) 


0.06% 


-0.08% 


-0.17% 


0.41% 


0.34% 


0.29% 


0.20% 


0.24% 


(50,500,200) 


0.63% 


-0.14% 


-0.42% 


2.50% 


2.22% 


1.97% 


1.47% 


1.71% 


(100,1000,200) 


0.47% 


-0.17% 


-0.48% 


2.12% 


1.88% 


1.62% 


1.20% 


1.41% 



TABLE III. 

Signal to noise ratio r in pp ^ tt + X at the LHC in the 2HDM with ^cp — 1- The 
combined results are for O;, O2 and O4. 



m^GeV) 


Ai 


T2 


/2 


Oi 




O2 


O4 


combined 


100 


-0.20% 


0.10% 


0.28% 


1.09% 


0.97% 


0.82% 


0.60% 


0.71% 


200 


-0.16% 


0.12% 


0.30% 


0.96% 


0.85% 


0.69% 


0.49% 


0.59% 


300 


-0.15% 


0.14% 


0.36% 


1.12% 


0.97% 


0.78% 


0.54% 


0.66% 


400 


-0.13% 


0.20% 


0.45% 


1.55% 


1.52% 


1.31% 


0.92% 


1.11% 


500 


-0.05% 


0.15% 


0.31% 


0.91% 


0.87% 


0.76% 


0.53% 


0.64% 
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TABLE IV. 

Signal to noise ratio r at the Tevatron and LHC with model independent top quark 
CEDM = 1 GcV^^ Qs = 1.97 x 10^^'^cm Qs, and the accuracies with which dj^ can 
be measured at the Tevatron and LHC with assumed numbers of events given in the 
text.df is given in unit of 10~^^ cm Qs- Oi is given only for leptonic events. The 
combined results are for OJ, O2, O4. 





T2 


/2 


Ol 


Ol 


01 


O2 


O4 


combined 


Tevatron 


r 


39.9 


40.9 


-35.2 


124 


88 


75 


47 


61 




10.9 


10.6 


12.3 


3.5 


4.9 


2.4 


3.0 


1.7 


LHC 


r 


-53.G 


-112 




238 


199 


151 


99 


128 




0.37 


0.18 




0.08 


0.10 


0.052 


0.064 


0.037 



TABLE V. 

The same as Table IV., but for d^. 







A2 


Ol 


o[ 


O2 


O4 


combined 


Tevatron 


r 


-80.5 


57.4 


214 


179 


169 


127 


146 


4 


5.4 


4.4" 


2.0 


2.4 


1.1 


1.1 


0.73 


LHC 


r 


83.2 


-21.1 


381 


332 


332 


248 


282 


4 


0.24 


5.8" 


0.052 


0.059 


0.025 


0.025 


0.017 



" Note that for ^2, only half of Nji can be used. 
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APPENDIX A: FORM FACTORS 



We give here the non-zero form factors for the matrix elements appeared in the 
text. They are written in terms of the conventional one-, two-, three- and four-point 



scalar loop integrals defined in Ref. . 

In the following, the form factors are given in the MSSM. 



i=i,2 



+F2{Co + Cii)(-p2, k, m-g, mt^,mt^)}^mg{-iy^H\cp 



(Al) 



jDB ^ J2 {Dis{-p2,PA,P3,mi mg,mq,mg)}^mg{-iy^^iXcp 
3=1,2 47r 



(A2) 



/1*1 = -2D^(P1-P4-P2-P4)/S 

ft = ^D^s , 



(A3) 
(A4) 
(A5) 



3=1,2 



ff = {Coi-P2,k,mg,mi^,mi.)}^mg{-iy^H\cp 



(A6) 



= E {-Boil ml ml) + B,{ml m% m? )} 



3=1,2 



X 



4:71 



m-gi-iy+HXcp/ii-m'',). 



(AT) 



vl,t 



E {-Coi-p2,PA, mi, , m-g, mg)}^^mg{-iy+^i\cp, 
3=1,2 '^^ 



Ar = 2D7(t-m^). 



(AJ 



ff'* = E {-Co{pi,-p3,mi^,mg,mg)}^^mg{-iy^^i\cp, 

3=1,2 



ff'' = ADy{i-m 



At'* = 2D7(t-m 



(A9) 

(AlO) 
(All) 



where 
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= E {FiCn(-P2,P4, mi.,m~g, rrig) (A12) 
+F2{Co + Cu){-p2,P4, rUg, mt.,mi.)}^mg{-iy+^iXcp 



j.bo.i,t ^ ^ {_4n2,}^mg{-iy+H\cpF, (A13) 

j^bo.l,t ^ ^ ^^2(^^ ^ 2L'i2 + ms + i^24 + i^25) " m^Do " 2^2 ' P4(2L'l2 (Al4) 

+i?24 - D26) - 2pi ■ p4(i^25 - D26) + 4D27}-^mg{-iy+HXcpFi 
/^''* = E {-4i?26}^m,(-l)^+^acpFi (A15) 
^60.1, ^ j2 {4p25 - i^26)}^m5(-l)^+^acpFi (A16) 

i=l,2 

i=i,2 47r 
/if = E {-4(^13 + i^26)}^m,(-iy+^zAcpFi (A18) 

/i3^''* = E {2D,2}^mg{-iy+HXcpF, (A19) 

j=l,2 

^60.1, ^ ^ {-2L>i3}^m5(-iy+^acpFi (A20) 

i=i,2 47r 

^.0.1, ^ ^ {-2Di3}^m^(-l)^+^zAcpFi (A21) 

j=i,2 47r 

/ir''* = E {2i^i2}^m,(-l)^+^zAcpFi (A22) 

The above D functions have the arguments {—p2-iP4.-iPz-i'n^ij-i''^9-i''^g-i^9)- 

fbo.2,t ^ ^ {_^D2r}^mg{-iy+HXcpF2 (A23) 

fbo.2,t ^ ^ ^ D2e)}^m~g(-iy+'iXcpF2 (A24) 

= ^ {4(^25 - i?26)}^m^(-l)^+^zAcpF2 (A25) 

i=l,2 

/9°"''* - E {4(^0 + 2Di2 + 1^24 - D2e)}^mg{-iy^HXcpF2 (A26) 
/lo"''* = E {-4pi3 + D2e)}^mg{-iy+HXcpF2 (A27) 

The above D functions have the arguments {—p2,P4,P3,mg,m£.,m£.,mf.). 
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fto.3,t ^ ^ {_^D,r}^m^i-iy^'tXcpF,F, (A28) 

^6o.3,* ^ ^ {4{-D,, + D,,)}^m-,{-iy^HXcpF,F2 (A29) 

^6ox3,* ^ ^ ^ /)^g)}^^,(_l).+i,A^pFiF2 (A30) 

.=1,2 47r 

^bo.3,t ^ j2 {A{D^, - + D,, - D2,)}^m-,{-iy^HXcpF,F, (A31) 
.=1,2 47r 

/lo"'' = E {-4(^11 - ^12 + /^2i - /^24)}^m^(-l)^+'acpFiF2 (A32) 
.=1,2 47r 

/is"''* = E {2(^12 - D,,)}^m-,{-iy^HXcpF,F2 (A33) 
.=1,2 47r 

f,o.3,t ^ ^ ^_2p^^ _ ^^^)}^ (_i).+i,^^^^^^^ (A34) 
.=1,2 47r 

The above D functions have the arguments {—p2,P4:, —Pi,nT'i ,''T^g,nT'g,nT'i )- 

3 1 ^ ^ 

In the above, the color factors Fi = - , F2 = — , and Xcp = 2Im{aib\) = 

2 6 

sin 29 1 sin Pt ■ 

The form factors of 2HDM can be obtained from those of MSSM by setting j = 
1 , Fi = 1 , F2 = , f^^ = /^^ = , and the following substitutions: 

Xcp ^ 2a~a = -2^cp , (A35) 
m-g^ rrit , nii^ m^^ , (A36) 



as \/2m^GF 



In 2HDM, the form factors of Fig.2 (n) defined in Eq.(^) are 



(A37) 



/r = ^^lf^ {im,[m',Co - P3 ■ p,{2C22 - 2C23 + Co) + ^]} , (A38) 
f2 = \/f^{^m,[Co + A{C22 - C23)]} , (A39) 



= ^^M££{4,„,c.} , (A40) 

where the C function arguments are (p4, —k,mt,mt,mt). 
The constant model-independent form factors are: 

D' = idt , (A41) 

= -2idt{pi ■P4-P2- Pi)/s , (A42) 

/s = ^^dt/s , (A43) 



21 



-fs /S 




J2 — ^^"t 5 


(' A /I 


f'"'i — AAA Iff ™2\ 






(A47) 




(A48) 


/i = -idt , 


(A49) 




(A50) 



where dt is defined in the text. 
Other definitions: 

k = Pi+P2=P3+Pi, s = k'^, i = q^ = {p2-pif, u={p2-pzf, (A51) 
= {-Pi + p^iYe^ ■ €4 + (2p4 + Ps) ■ £364 - (2P3 + P4) ■ e^e"^ . 



APPENDIX B: SOLUTION OF THE NEUTRINO MOMENTA 

We give here briefiy the method of solving the neutrino momenta in Eqs.(^5]). Let 

Pi+ = Pu Pi- = P'l, Pb = P'b, Pui = {E,X,Y,Z), p^^ = {E',X',Y',Z'), pib = pi +Pb, 
P'lb = P'i+ P'h-, Po = -{Pb + p'b + Pi+ P'l)- Then we can get 

X = aE + hZ + 5 , Y = cE + dZ + i , (Bl) 
X' = a'E' + b'Z' + 6' , Y' = c'E' + d'Z' + ^ , 



where 



EibP^ - EipI -pffcpf + p!pI 

a = , h = , (B2) 

_ -EibPf + Eipl pff^pf - pIpI 

-{mj - pl)/2p\ + ml,/2pl {mj - pi) / 2pf - ml, / 2pl 

' = A^ ' ^ = A^ ' 

Ai = pfbP^ - plpt , 

and a',b',c',d',6',^' are obtained from the above equations by the substitutions of 
momenta without ' to those with '. We further express E, E' in terms of Z, Z': 

E = fZ + gZ' + h , E' = fZ + g'Z' + h' , (B3) 

where 
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-bc' + da' -h'c' + d'a' , {^p^^-5-5')d + a'{pl-i-^') 

^ = A, ' ^ = ' ^^'^ 

, _ -ad + cb I _ -ad' + ch' , _ {pi - i - i')a + c(pg - 5 - 5') 
^ " A2 ' ^" A^ ' A^ ' 

A2 = ac' — ca' . 

Inserting the above expressions into E"^ = + + and E'"^ = X'^ + F'^ + Z'^, 
we get the following two quadratic equations for Z, Z': 

AZ^ + BZZ' + CZ'^ + D + UZ + VZ' = , (B5) 
A'Z'^ + B'ZZ' + C'Z'^ + + f/'Z' + = , (B6) 

with 

A = /2 _ 1 _ (a/ + 6)2 - (c/ ^d)\ B = 2[fg - ag{af + b) - cg{cf + d)] , (B7) 
C = g^- a^g^ - c^g^ , D = - {ah + 6)^ - {ch + 0' , 
U = 2[fh - {ah + 6){af + b)- {ch + ^){cf + d)] , 
V = 2[gh - ag{ah + 6) - cg{ch + ^)] , 

and A', B', C, D', U', V can be gotten from the above equations by replacing the vari- 
ables without ' with the ones with ' and interchanging /', g'. 

Eqs.([B5D (PB|) can obviously give at most four real solutions for {Z,Z'). They can 
be solved by standard methods: solving Z' in ( |B5| ) as functions of Z, we get two 
expressions for Z', then inserting them separately into ( p6|) , we shall get the same 
quartic equation of Z which has up to four real solutions. Although those solutions 
will be doubled when we calculate Z' with the two expressions obtained from 
half of them are false solutions to fB6f) . 



APPENDIX C: ANALYSIS ON THE CONTAMINATIONS TO CP 
VIOLATING EFFECTS IN PP COLLISION 

We assume the partons inside the proton have no intrinsic transverse momenta and 
denote the parton level reaction as a{xi)d{x2) —>■ ti with xi and X2 being the Bjorken 
scale parameters. The hadronic level cross section can be written as 

da = 2fP{xi)fi{x2)da{ad tt)dxidx2/{l + S^-a) (CI) 
= W{^l)mX2) + fl{x2)fi{x,)] + m{x,)fi{x2) - mX2)fi{x,)]} 
xdcr{ad ^ tt)/{l + Sa,a)dxidx2, 
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where f^, /| are the parton distribution functions of a and a in proton. The function 
X2) — fa{xi)fi{x2) -\- fl{x2) fi{xi) is CP-even(its CP transformation is just the 
interchanging of Xi, ^2), so that no contamination comes from it. Now we look at 
the function F~{xi,X2) — fa{^i)fi{^2) — /a (^2) /I (2^1)- Contaminations should come 
from F~{xi,X2) term since it is CP-odd. F'~{xi,X2) — Q lor a — g. There are no 
contaminations from initial gluon reaction. All contaminations come from the initial 
quark reactions which are subdominant processes at the LHC. There are non-zero 
CP violating contaminations only when the observables contain asymmetry between 
xi and X2- Since the Lorentz invariant observables can always be calculated in the 
center of mass frame of the parton which depending only on ^1X2, they will receive 
no contaminations. Therfore, within the framework of parton model, there are no 
contaminations to Lorentz invariant observables. In our study, they are /2,0i,0;, 02,04. 
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Figure Captions 



Fig.l Feynman diagrams of tree level and 1-loop SUSY QCD corrections of CP violation 
to qq tt 

Fig. 2 Feynman diagrams of tree level and 1-loop SUSY QCD corrections of CP violation 
to gg tt 

Fig. 3 Feynman diagrams of top quark CEDM corrections of CP violation to qq — > tt, 
gg tt 
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